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Abstract 

The new version of the parton shower module APACIC++ for the SHERPA event 
generator framework is presented. It incorporates some features, that are specific 
for the consistent merging with multi-particle matrix elements at tree-level. This 
publication also includes some exemplary results and a short description of the 
upgraded class structure of APACIC++, version 2.0. 
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1 Introduction 


Event generators for the simulation of high-energetic reactions of particles play an im¬ 
portant role in understanding inclusive and exclusive final states at collider experiments. 
In order to live up to the complicated task of a complete simulation of such processes, 
event generators decompose the full reaction into a sequence of individual steps, usu¬ 
ally divided by a characteristic scale. Typically, the starting point for the simulation 
of an individual event is the selection of the hardest process in it, commonly this is 
denoted as the signal process and described at the parton level. The emerging partons 
are then passed to the parton shower. This part of the simulation describes additional 
multiple parton emission. In modelling this Bremsstrahlung-like radiation, certain ap¬ 
proximations are made, that retain the leading soft or collinear contributions and allow 
the description of complex radiation patterns through individual parton splittings at 
decreasing scales. When the splitting scale reaches a cut-off value, typically of the order 
of a few Aqcd, the parton shower stops. At these scales, the realm of perturbative 
QCD is left, and the regime of soft QCD is reached. There, the final translation of the 
parton to primordial hadrons is modelled through phenomenological approaches. If the 
primordial hadrons are unstable, they are further decayed until only stable hadrons are 
left. 

In such a framework, the parton shower provides the link between perturbatively cal¬ 
culable differential cross sections at the parton level and models for their transition to 
observable hadrons with phenomenological parameters, which need to be tuned to data. 
The inclusion of the parton shower does not only create more realistic high-multiplicity 
parton final states from low parton multiplicities available through corresponding exact 
matrix elements, it also reduces the average distance of the partons in momentum space 
down to a fixed size and thus ensures that the hadronisation parameter tunes are rather 
independent of the hard process in question. It is this aspect of the parton shower that 
renders it indispensable for a meaningful simulation with any predictive power. This can 
then be used in order to calculate hadronisation corrections or to judge detector effects 
in certain processes. In most cases, results from parton shower Monte Carlos have been 
in astonishing agreement with data, giving rise to the conhdence that simulation tools 
can also be used to predict signals and backgrounds for the current and the next round 
of collider experiments. 

Although, for more than two decades, parton showers have been widely used in multi¬ 
purpose event generators, details of their specific implementation traditionally depend 
on certain choices and, due to recent developments, on the specific form the hard matrix 
elements generate the initial parton configurations. Traditional multi-purpose event 
generators, like Pythia nia or Herwig PH usually have a 2 ^ 2 process at leading 
order (LO) as the signal process. There, after dehning the corresponding, process- 
dependent starting conditions for the shower, the latter is allowed to evolve freely. Only 
recently, this treatment has been consistently extended to next-to leading order (NLO) 
precision for 2^2 processes in the MCONLO framework iniinj, which uses the machinery 
of Herwig. There, the NLO calculation is modified in order to match the requirements 
imposed by the specihc form of the parton shower. Another approach, going beyond 
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2 —i> 2 processes, has been taken in SHERPA [7j; there, a fully automated consistent 
merging of 2 —>■ n processes at the tree-level with the parton shower according to the 
formalism of IHlIi has been implemented and tested P1‘. For the merging with such 
multi-leg tree-level matrix elements, the parton shower has to be supplemented with the 
determination of more involved starting conditions and with additional constraints on 
the phase space of the emitted particles. 

In this paper the specihc realisation and implementation of the parton shower in SHERPA 
will be presented. A first version of this module, APACIC++, version 1.0, has been 
published some time ago in h covered the parton shower in the hnal state only, in¬ 
cluding some algorithms for the merging with the matrix elements provided by AMEGIC++ 
[HI EH). In this hrst version, many of the steering and service classes were still located 
within APACIC++; as time went by, APACIC++ transformed from a stand-alone code to 
a mere module of the full framework SHERPA. In its present state, version 2.0, which 
is discussed in the following, APACIC++ thus includes the parton shower in the initial 
and hnal state and an improved handling of algorithms necessary for the parton shower 
aspects of the merging procedure^. The outline of this paper is as follows: After briehy 
reviewing the parton shower formalism in Sec. |2l the ideas underlying its particular 
implementation in APACIC++ will be discussed in Sec. El There, special focus will be 
on those aspects, which are specihc for the merging procedure with matrix elements. 
In Sec. m exemplary results will be presented, which compare the parton shower results 
of APACIC++ with analytical approaches and with data. Details of the implementation, 
including a short description of the class structure of APACIC++ will be given in the 
appendix, App. O 

2 Parton shower 

This section comprises a brief overview of the parton shower formalism and the corre¬ 
sponding algorithms implemented in APACIC++. There, the sequence of parton emissions 
in the shower evolution is organised by virtuality as the ordering parameter. In this re¬ 
spect, the algorithm presented here is closely related to the implementation in Pythia 
El 12]^. However, the APACIC++ version differs in some details like the treatment of 
massive particles and particular scale choices for the evaluation of coupling constants 
and parton density functions (PDFs). In addition, it incorporates some unique features 
that facilitate the merging of matrix elements with the parton shower according to the 

similar approach for the merging of matrix elements with the dipole cascade in has been taken 
in m; for some specific processes the formalism has been used to merge multi-jet matrix elements with 
Pythia and Herwig m A slightly different method m has been advocated by the authors of Alpgen 
who merge their matrix elements with the parton shower of Herwig. 

^However, it should be noted that situations like deep-inelastic scattering are still beyond the reach 
of APACIC++ and, thus, of SHERPA. 

^Note that recently, a reformulation of Pythias parton shower has been presented m, which employs 
transverse momentum as ordering parameter. It has been implemented into the recent version of Pythia 
|2()|. In 123 an alternative evolution variable, also related to transverse momentum, has been discussed. 
It has already been implemented into Herwig++ \22\ . 
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formalism of iHiinj- 


Basics of part on showering 

The parton shower evolution relies on the fact that parton emission processes become 
singular in the soft or collinear limit. When the available phase space is cut accordingly, 
these singularities translate into large logarithms, which can be resummed according to 
the DGLAP evolution [23 12111211 1^ - By taking into account the leading logarithms 
only, the parton shower picture reduces complex radiation patterns of multiple parton 
emissions to chains of individual independent parton splittings. They are organised 
in a probabilistic manner by an ordering parameter, usually some quantity like the 
virtual mass of the decaying parton or the transverse momentum of the decay products. 
Suitable constraints on this ordering parameter avoid singular regions of phase space. 
The probability for no parton splitting to occur between two scales to < G is encoded 
in the Sudakov form factor. It is given by 

Ta—>bc('^) / ) 


Aa(to,ti) = exp 


dt 

t 


dz 


27r 


^0 


E 

6,c 


where p± is the transverse momentum. It can be written as a function of the scale t and 
the energy splitting variable z. Pa-,bc{z) denotes the splitting function for the branching 
a ^ be. A complete list of splitting functions implemented in APACIC++ can be found 
in App. IHl Because of its interpretation as no-splitting probability between two scales, 
ratios of two Sudakov form factors like 


P 


nobranch 


A(to, ti) 

A(toG) 


( 2 ) 


yield the probability for no emission between ti and t, which could be resolved at the 
scale to- Given a hard starting scale ti, this allows the scale for the actual branching t 
to be generated by equating a random number R with this ratio and solving for t. With 
a second random number then the splitting variable z can be selected. It is distributed 
according to 

Upn, = PU ). (3) 

This probabilistic interpretation allows the formulation of the parton shower as a Markov- 
chain of independent 1 —>■ 2-branchings, where the scale ti, at which the former splitting 
occurs, sets the upper limit for the subsequent branching. Organising the parton shower 
in such a way results in what is known as “forward evolution”. This type of evolution 
is employed for parton showers in the final state, i.e. for the “time-like parton shower”, 
where the resulting parton ensemble is not subject of any other constraint. 

For the initial state, however, i.e. for space-like parton evolution, the situation changes 
drastically. This is because standard forward evolution of a parton ensemble from some 
comparably low fixed hadronic scale Qq < 0, of the order of a few Aqcd, to the fixed scale 
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of hard interaction, Ql, distributed according to the appropriate matrix elements, 

would be highly inefficient. Therefore it is more convenient to start with specifying the 
scale of the core process and its momentum fractions xi and X 2 , and subsequently 
evolve “backwards” to the partons assumed to be resolved from the incoming hadrons. 
The correct way to perform this backward evolution has been introduced in j2Zll2Hl; it 
boils down to modifying the Sudakov form factors by dividing by the appropriate PDF at 
the given scale and at the relevant Bjorken-a:. The resulting no-branching probabilities 
are then given by^ 


, _ fix, t) 

nobranch — n / , ^ 

fix, ti, 


^jto, ti) 


(4) 


Again, equating this probability with a random number and solving for t dehnes the 
previous branching scale. This time, however, the corresponding splitting variable ^ is 
distributed according to 


-Psplit — 


asip±) Pjz) ^ {X \ 
2ti z •' U ’ ) ’ 


(5) 


where x is the momentum fraction of the decaying parton, and z = x/xi with xi being 
the momentum fraction of the resulting parton from which the backward step started. 
In both time-like and space-like evolution, the notion of soft colour coherence |281ld()ll!?n 
(23 EH plays a crucial role. It results in an angular ordering of subsequent branchings, 
which can be implemented directly through a corresponding choice of the evolution 
variable. This has been done in Herwig 13111, where a suitable angular variable has 
been chosen. Alternatively, colour coherence can be implemented by choosing transverse 
momentum as ordering parameter; this choice has been made in Ariadne m which bases 
its multiple emission treatment on splitting colour dipoles. Recently, two new formalisms 
to incorporate ordering according to transverse momentum into the more conventional 
parton picture have been presented in US HH- In its most trivial version, however, 
angular ordering can be imposed through a direct veto on increasing opening angles. 
This is how it has been implemented in Pythia, where virtuality is the evolution variable. 
Of course, such a way of implementing colour coherence is by far not as sophisticated 
as the methods above; nevertheless this method has also been chosen in the framework 
of APACIC++. 


Variables in APACIC++ 

As indicated above, the parton shower in APACIC++ is organised in terms of virtuality, 
i.e. in terms of the virtual mass of the decaying particle, and the splitting variable 2 ; is 
interpreted as the energy fraction of one decay product (daughter) w.r.t. the decaying 
parton (mother). 


more efficient, non-Markovian, evolution algorithm is currently under investigation m- 
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1. Time-like evolution: 


For parton branchings a ^ be in the hnal state, the evolution variable t is given 
by 

t = ta=pl = {Pb+Pcf ■ (6) 

The splitting variable 2 ; is dehned as energy fraction, i.e. 


taken in the rest frame of the complete hnal state parton shower^. For the respec¬ 
tive transverse momentum, which serves as an argument in the running coupling 
constant and as a low-energy cut-off, a dehnition following has been chosen. 
There, 

kl = 2Min{E2, e2}(1 - cose ,,), (8) 

which in terms of ta and z translates into 

= Min (9) 


There is, however, some residual freedom in the exact choice of scales, which can 
be altered by a (logarithmically) small factor. In APACIC++, this freedom is used 
to define the transverse momentum as 

( 10 ) 


Demanding that all p± are larger than some minimal cut-off, pf\ immediately 
poses a constraint on the allowed (ta, 2 ;)-range for individual branchings, namely 


ta> to = 


and 2 ; G 


^0 

ta + to 


ta + to 


( 11 ) 


In addition, the transverse momentum must be positive. Using the dehnition of 
the splitting variable, Eq. dZD, yields a kinematical condition on z, namely 


2; e 


1 - 


I tg 
E2 


1 
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1 + 


hi. 


( 12 ) 


For the explicit angular veto employed to model coherence ehects, the following 
approximation for the angle is used 


'be 




(13) 


^In APACIC++, the time-like evolution of a given parton ensemble is performed in its rest frame. 
After the shower terminates, the resulting partonic final state is fully reconstructed and boosted into 
the relevant (lab) frame. 
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It stems from ta = 2z{l — z)E‘^ (1 — cos 6bc)- This approximation is valid for small 
9bc, but it can conveniently be used for large angle emission too, since the actual 
value is irrelevant for its role as ordering parameter. 

Having split parton a into b and c, the two offsprings are massless hrst. They gain 
a virtual mass by continuing the parton shower evolution for another step with 
starting scale ta- Of course, the sum of their virtual masses cannot be larger than 
the virtual mass of a, i.e. 

■y/ffe + \/^ < \/^ • (14) 

This, together with the angular ordering constraint, sets some additional limits 
on the splitting variables. However, there is still an issue to be resolved: When 
^ was defined, it was implicitly assumed that both outgoing partons are on their 
mass-shell; this is not true any longer, which implies that the branching kinematics 
of a —>■ be has to be redehned. In APACIC++, this is achieved through a redehnition 
of the energy splitting variable 2 ;, also quite along the lines of what happens inside 
Pythia. In both codes, the modihed splitting variable reads 

~ f 1 \ tc) ta tb tc 

where 

A(a, b,c) = ^/(a — 6 — c)2 — 46c. (16) 

This modiheation constitutes a simple linear transformation from the massless to 
the massive domain, which is always possible as long as the condition given by Eq. 
m is satisfied. 

If one or more of the partons a, b, and c are massive, the dehnitions above change 
as follows: 



• The evolution variable is now dehned as ta = ta — This in addition to 
the usage of mass terms in the slitting kernels (cf. App.^, leads to a slightly 
modihed Sudakov form factor. 

• The argument of the running coupling is now given by the “complete” k± 
formula, i.e. it is dehned through 


P± = 




(1 - cosObc ), 


which assumes on-shell, but not necessarily massless daughters. 


(17) 


The other conditions introduced above are, in principle, not altered, provided that 
the kinematical variable ta is used rather than the evolution variable ta- Conse¬ 
quently, the kinematical z-domain is still given by Eq. m, and z is translated 
into the massive case by Eq. (USD. This leads to a reinterpretation of the splitting 
variable, which is necessary, since the splitting functions are calculated with z 
being the light-cone momentum fraction. 
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2. Space-like evolution: 


For the parton shower in the initial state, i.e. space-like showers, the evolution 
proceeds backwards. The ordering parameter here is the virtual mass of the re¬ 
spective initial state parton. Therefore in branchings b ^ ac the scale is given by 


t = 'ta=pl<0, ( 18 ) 

and, wherever it is needed, the absolute value is taken. The dehnition of the 
splitting variable is a little bit more tricky. In order to ensure four-momentum 
conservation in the reconstruction of the showering kinematics, the two shower 
branches of the two incoming partons are conpled. For the two incoming particles, 
labelled by a and a, the Sndakov form factors, cf. Eq. dH), are used to choose their 
virtual mass, ta and ta- The parton with larger off-shellness is selected for the 
reconstrnction of the corresponding backward step. Assume that ta < ta, and, 
hence, that a is selected. Then z is defined throngh 


{Pa+Pdf 
{Pb + Pd)"^ ’ 


( 19 ) 


clearly a Lorentz-invariant measnre. This implies that, step by step, the c.m. 
energy sqnared in the parton system is enhanced by a factor \jz. In the same 
way, the Bjorken-x of partons a and b are related by 


Xb = xjz , 


( 20 ) 


where the hrst Bjorken-x of the partons entering the hard process are taken directly 
from the matrix element evaluation. In this respect, at each step of the space-like 
shower evolution, the current pair of partons is oriented along the beam axis, with 
Bjorken-x as if they were massless. It is this Bjorken-x and a scale = p\ that 
constitute the parameters for the calculation of the PDF at the corresponding 
step. In the c.m. system of partons b and a, the transverse momentnm of a and c 
is given by 

kl = -{l-z)ta. (21) 

Following a similar reasoning as in the hnal state treatment, the transverse mo¬ 
mentnm entering the rnnning conpling as well as the PDF evalnation is defined 
throngh 

Pi = -^-^ta- ( 22 ) 

Having thus constructed the backward step b ^ ac leading to the parton b, it is 
clear that c may experience a final state shower evolntion. Its starting scale, by 
default, is given by ta- For kinematic reasons, however, the actual scale that is 
chosen in the shower evolution must also satisfy 

^ _ ^ad^bd Sad^bd ^ ^ ^ 

^ ^c,max r “T 
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with 


^ad ^a.d ta ta ; 

^ad ^{^adi ^0,1 


s'hd = Sbd -tb-ta ( 24 ) 

\d = Ks'bd^h,t~a) ( 25 ) 


where Sad = {Va + 'PdY and Sbd = {Vb + Pa)^ = Sad/z. Similar to the approximation 
of Eq. (Uni), the opening angle of this splitting can be estimated by 


e 


ac 



( 26 ) 


Knowing tb and tc as well as the four-momenta pa and pa, the four-momenta pb 
and Pc are constructed explicitely in the c.m. frame of pa and pa- In this system, 
energies and momenta are fixed according to 


Ec 


^bd 4“ 4“ tc 



( 27 ) 


Pz,c 

P±c 


^bd ^E^Ec 



( 28 ) 

( 29 ) 


where the transverse momentum is distributed uniformly in its azimuthal angle. 
Furthermore, 


Pb=Pa+Pc- ( 30 ) 

Subsequently, the system is boosted into the rest frame of partons b and d and 

rotated such that their momenta point along the beam axis. 

Colour treatment 

When the parton shower terminates, the resulting multi-parton ensemble needs to 
be transformed into hadrons. Due to a lack of quantitative understanding of non- 
perturbative physics, this is achieved through phenomenological models. These models 
have some different underlying physics assumptions, highlighted by the examples of the 
string model ISni EEl E3 EHl EHl EDI stressing the role of colour coherence EH, and of 
cluster fragmentation models [iniEnilI3llI3E3 that are closer to the independent frag¬ 
mentation approach EHl and the concept of pre-confinement m- However, they have in 
common that they rely on having as input a parton ensemble with a well-dehned colour 
structure in the Nc ^ oo limit. This necessitates that the parton shower must distribute 
colours in this limit. In most cases, like, e.g. in a splitting q ^ qg the colour structure 
is unambiguously defined, ensuring that a well-defined colour structure at the beginning 
of the parton shower evolution can be mapped onto a well-defined colour structure at its 
end. However, there are cases with ambiguities, namely in splittings of the type g gg. 
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In this case, there are two ways of locally distributing colour. Denoting a colour state 
by its colour triplet and anti-triplet labels in the Nc ^ oo limit, (a, b), these two ways 
in the splitting 5 f(l) —> g{2)g{3) can be visualised as 

(a, b) —> (a, c) -f (c, b) or (a, b) —>• (c, b) + (a, c). (31) 

In APACIC++, the choice is made in the following way: In the shower evolution, the 
decaying gluon 1 was produced in another parton branching process, where the other 
decay product is colour-connected to the gluon. For both gluons 2 and 3, the transverse 
momentum w.r.t. this colour partner is determined according to Eq. ©. The choice is 
then made such that the gluon with the smaller relative k± is colour connected with this 
partner parton. 

Initialisation of the parton shower 

Within the SHERPA framework, an algorithm along the lines of IHlIi for the consis¬ 
tent merging of tree-level matrix elements for multi-particle production with the parton 
shower of APACIC++ has been implemented |1H]- The key idea of this algorithm is to 
separate the phase space for parton emission into the hard region of jet production ac¬ 
counted for by suitable tree-level matrix elements and the softer region of jet evolution 
covered by the parton shower. Then, extra weights are applied on the former and ve¬ 
toes on the latter, such that the overall dependence on the separation cut is minimal. 
The separation is achieved through a k± measure IMl US EO]- The weight attached 
to the matrix elements is constructed employing Sudakov form factors, thereby taking 
into account those terms that would appear in a corresponding parton shower evolution. 
Therefore, a “pseudo parton shower history” is reconstructed by clustering the initial 
and hnal state particles from the tree-level matrix element according to the k± algo¬ 
rithm. This procedure provides the scales necessary for the evaluation of the weight, 
namely the nodal values of the different k±, where two jets have been merged into one. 
Within APACIC++, the reconstructed “shower history” is utilised for the determination 
of initial conditions for the shower evolution. Then, starting from the hard 2 —>■ 2 
core process, all partons obtain as starting scale for their shower evolution the nodal 
virtuality, where they emerge for the hrst time. It should be stressed at this point 
that there is some residual mismatch in the parton shower variables used in APACIC++, 
namely virtuality t = and the scales {k±_ measures) used in the Sudakov form factors 
employed in the reweighting of the matrix elements. This point will be discussed in 
more detail in the following section. 

For the scale of the four partons entering or leaving the 2 —> 2 process, there are a few 
options, dehning the respective choice of scale fstart^ 

• Processes of the type > qq\ fstart = s = {pq + Pq)"^ ■ 

• Drell-Yan type (e.g. qq e+e"): tstart = s = {pq+ PqY- 

• “Drell-Yan-|-jet”-type (e.g. qq W^*'^g)\ For the incoming partons the choice is 
4tart = s = {pq + PqY- For the hnal state parton, however, the start scale for the 
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shower evolution is different, namely the nodal value where it has been produced, 
usually t or it. 

• QCD processes: In this case, different colour flows are usually competing. The 
winner is selected according to the respective contribution, which is related to a 
specihc propagator structure (s-, t-, or u-channel). The choice of starting scale 
follows the selection of the colour structure; in most cases therefore the starting 
scale tstart is the minimum of t and it. 

In order to account for colour coherence, a maximal allowed branching angle dcr\t has 
to be determined for each parton. Using the knowledge of the colour connections of 
the 2 —2 process, Oait is dehned as the angle between the parton in question and the 
parton colour connected to it. In case of gluons the choice of this angle is, of course, 
ambiguous; there, the maximum of both possible values is taken. 

Having fixed the starting scales of the hardest partons, the starting scales for the evolu¬ 
tion of softer partons are easy to obtain. In principle, there are only two ways, in which 
softer partons may emerge, which lead to a slightly different treatment. 

• Branchings in the hnal state: 

There the harder, i.e. the more energetic, of the two offsprings inherits as starting 
scale and angle the values of the decaying parton, whereas the parameters of the 
softer offspring are taken directly from the node. In other words, the initial virtual 
mass equals the invariant mass of the pair, and the starting angle is given by its 
opening angle. 

• Branchings in the initial state: 

There, in analogy to the treatment of the space-like shower evolution, the virtual 
mass and the starting angle of the time-like offspring are given by the virtual mass 
of the initial state parton that branches “backwards” and the respective opening 
angle w.r.t. the corresponding beam. The “new” initial state parton inherits the 
starting scale and angle of the line pointing towards the “core” process. 

Merging issues 

Having fixed the starting scales, the parton shower can be evolved. Due to the merging 
prescription, emissions inside the parton shower have to be vetoed, if they result in the 
production of partons with transverse momenta k\_ larger than a given jet resolution 
scale, p™*. Parton showers attached to matrix elements with the highest multiplicity 
must live up for the production of eventual extra jets. In this case, the veto scale is 
not the jet resolution scale but rather given by the smallest k± of the partons stemming 
from the matrix element. 

However, in the parton evolution of such multi-parton hnal states, there is a last subtlety. 
It is connected to the fact that, through showering, the partons may acquire a virtual 
mass different from their on-shell mass used in the matrix element. This recoil has to 
be compensated for, if possible in such a way that the other partons stemming from the 
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matrix element are least perturbed. This is solved in the following way: The merging 
procedure results in the reconstruction of a pseudo parton shower history of the parton 
ensemble. Thus, to each parton produced in the matrix element a partner can be 
assigned such that both emerge in a branching of the pseudo parton shower history. If 
one of them or both acquire a mass, the respective partner takes care of the recoil. The 
strategy employed for this is very similar to the one used when the initially massless 
partons in the showering acquire a mass; in both cases the energy splitting variable is 
shifted through a linear transformation, for the shower case, cf. Eq. (USD. In the initial 
branchings of the shower initiators, there are two cases to be considered: 

• Both partners experience a parton shower evolution. Then, the reconstructed 
energy splitting variable of the branching where they have been produced, is shifted 
according to Eq. (USD- 

• In the branch a —> be, only one partner (b) experiences a parton shower evolution, 
the mass of its partner is (c) fixed, since it can be thought of as an internal line of 
the matrix element - it branches according to the pseudo parton shower history. 
Then z is shifted according to 

+ X{ta,h,tc) ta+tb-tc 

where denote the matrix element masses and ti the new virtualities. This 
modification results in a change of the opening angle 9bc, leading to a modified 
momentum even for the particle with the mass unchanged. In order to comply 
with four-momentum conservation, its offsprings have to be adjusted, too. This 
transformation can easily be achieved through a rotation followed by a boost along 
the direction of particle c. 

Within APACIC++, an additional veto on “losing” a jet produced in the matrix element 
is introduced. Losing a jet through the shower is possible due to a mismatch of the 
quantities responsible for jet definition {k±) and jet evolution (virtual mass). This 
leads to recoils imposed by the shower evolution resulting in a change of k±, while the 
virtuality of the internal line in question is preserved. Therefore, a check on the number 
of jets is performed after the shower evolution is hnished. This is done by demanding 
that the - now off-shell - shower seeds are still separated in k±. If this is not the case 
the event is rejected, and a new event is generated from the beginning, i.e. with a new 
kinematical situation but the same flavour constellation of the jets. 
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branching 


Knot 



Figure 1: The basic building blocks of shower emission are binary branchings, repre¬ 
sented by Knots. 

3 Implementation 

This section focuses on the basic strategies according to which the formalism discussed 
in the previous section, Sec. El is implemented in APACIC++. A more detailed reference 
to all individual classes can be found in App. El 

The basic unit for the realisation of the parton shower in terms of a computer program is 
the individual branching of a single parton. Connecting such branchings automatically 
leads to the Markov structure of the full emission pattern. In other words, the task of 
a parton shower program is to hll and to connect such individual branchings. It is the 
strength of an object-oriented programming language like C++ that it is well suited to 
map the underlying physical pictures onto program code by using an appropriate class 
structure. 

Representation of the parton shower 

Following the reasoning above, the basic structure representing a single branching is a 
Knot, reflecting the a ^ be binary decay structure inherent to the parton shower. A 
Knot carries information on the incoming particle a, encoded in the class Particle, and 
on the Knot, where it originates from. When it is hlled, i.e. when the decay is specihed, 
it also yields the two offsprings and their respective decay Knots. In addition, a number 
of other quantities are stored in a Knot, namely 

• the scale 

• the energy splitting parameter z, 

• the squared energy, 

• and its minimal virtuality t^a\ 

• the respective opening angle in the approximation of Eqs. (naEni) ^bci 

• the azimuthal angle cj), 

• and the Bjorken-x of the particle Xa (if needed). 
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Figure 2: Sketch of the mapping between radiation processes and the corresponding 
classes. The full radiation pattern is identihed as a chain of 1 ^ 2 processes, a Markov 
chain, which translates into the class Tree. A Tree in turn is realised as a list of linked 
Knots. The shower evolution of an event is represented by three Trees, one for the hnal 
state shower, and two for the initial state shower. 


These Knots are then linked in terms of a (binary) Markov chain, represented in the 
class Tree. Each Tree contains a pointer to its hrst Knot, the root Knot®. Starting 
from this root Knot, all other Knots are then accessible by successively following the 
pointer inside the Knots, thereby spanning the full Tree structure. Each event consist 
of three Trees, as indicated in Fig. El one Tree gives rise to the complete hnal state 
shower, eventually with a dummy particle as root Knot. The other two Trees represent 
the initial state showers from both sides, therefore they have one of the two particles 
entering the hardest subprocess as their respective root Knot. Since the full parton 
shower evolution is represented by these Trees, they provide routines to create and 
delete individual Knots, or to boost and rotate the full structure. The latter option is 
relevant for the transparent implementation of the initial state parton shower. 

Filling the Trees 

The Trees are now to be hlled by the corresponding parton shower. In APACIC++, two 
classes are responsible for this, namely the Final_State_Shower and the Initial_- 
State_Shower, who have access to one or two Trees, respectively. Since the latter, 
the Initial_State_Shower, produces particles that may undergo a time-like shower 
evolution, it also has a pointer to the former, the Final_State_Shower. Apart from 
switches steering the specihc way the shower works however, both showers need infor¬ 
mation from Sudakov form factors and they have to construct the respective time-like 
or space-like kinematics. These necessities are encoded in corresponding classes, named 
Timelike_Sudakov or Spacelike_Sudakov, and Timelike_Kinematics or Spacelike_- 

® Since a parton shower history is reconstructed in the merging procedure, this root Knot may reflect 
the already fully defined splitting of an internal line of some matrix element. 

^There are some options for specific questions made available in APACIC++, that go beyond the 
standard settings discussed above. 


14 









owns and employs object -employs object 


Figure 3: Relation of the main classes of APACIC++. 


Kinematics. It is clear that the Spacelike_Kinematics must use methods from its 
time-like counterpart, hence it has a pointer to Timelike_Kinematics. Since the two 
Kinematics-classes are responsible for jet vetoes, they have a pointer to a Jet_Finder 
class made available through the SHERPA framework and some flags steering its proper 
usage. An overview over the basic relations between the main classes of APACIC++ is 
depicted in Fig. 121 

In contrast to the Kinematics classes, the two Sudakov classes are a bit more intricate. 
Since most of their actions are related to the selection of splitting functions and their 
usage, both are derived from the class Splitting_Group, and contain a Splitting_- 
Group for each flavour, the shower can handle. Then, each Splitting_Group builds a 
wrapper around all splitting functions Pa^bd^) for all branchings allowed for a certain 
incoming particle a. As an example, consider the case of a gluon. It may undergo either 
9 ^ 99 OT & g ^ qq branching, each of which is represented by its corresponding 
Splitting_Function. Consequently, the class Splitting_Group responsible for the 
gluon branching contains one Splitting_Function for the gluon final state and one for 
each quark flavour. Any individual Splitting_Function incorporates information like 
its incoming and outgoing flavours, an estimate for its integral over z, and methods to 
extract z distributed according to Packed). In addition, Splitting_Groups allow to 
select a splitting mode and the corresponding flavours. 

The two Sudakov classes construct all such Splitting_Groups for physical branchings. 
When a parton, encoded in Knot is to be split, the respective Sudakov form factor 
implementation selects the physically relevant Splitting_Group, i.e. the one with the 
flavour of the parton as incoming or outgoing particle, and uses its estimated z integral 
to select a corresponding t Then vetos are applied, either in the procedure of the 
hit-or-miss method used in the Sudakov classes, or in the framework of the merging 
procedure to reject unwanted jets to be produced in the parton shower. The latter test 
is performed with the help of the corresponding Kinematics class. It should be noted 
that the Spacelike_Sudakov classes contain a PDF due to the backward evolution 
described above; in view of the possibility to have two different beam particles with 

® Other information needed for this selection, like the start scale fstart and the maximally allowed 
angle Aut are obtained from the Knot. 
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different PDFs the class Initial_State_Shower has two of them. 


4 Results 

In this section some results will be presented that validate the implementation of the 
parton shower algorithms, discussed in the previous sections. 

Comparison with analytic Sudakov form factors 

First of all, the parton shower of APACIC++ will be confronted with known analytical 
results dealing with the resummation of large logarithms emerging in multiple parton 
emission. Such results are available mainly for e’''e“ annihilation into hadrons. There, 
analytical calculations describe the relative rates for different jet multiplicities in the k±_ 
algorithm. For the massless case, results can be found in [^, massive quarks are dealt 
with in inn. All results presented here relate to e’''e annihilations at i?c.m. = 91.2 GeV, 
i.e. to the LEP 1 energy. 

To gain some hrst idea about the behaviour of the merged cross sections and the effect 
of the Sudakov form factors, consider Fig. |3J There, tree-level cross sections for the 
production of up to four jets, out of which two may be 6-jets (left) are contrasted with the 
weight applied to them in the merging procedure (right). The jets are defined through a 
Durham jet measure with resolution parameter y, which represents the minimal relative 
transverse momentum of two separated jets. Already here, it becomes apparent that 
the massive 6-quarks are suppressed w.r.t. massless ones in the cross section, whereas 
the weights are nearly identical for massless and massive jets. In Fig. Elthe merging has 
been performed; in the left plot the resulting jet cross sections and their sum are shown 
for the massless case, the right plot exhibits the corresponding cross sections when two 
6 quarks are involved. It is surprising how stable the total cross sections are w.r.t. the 
jet definition cut. In Fig. Elthe resulting two- and three-jet rates are displayed. Clearly, 
due to the “dead cone” effect the 6-quarks tend to radiate less, leading to an enhanced 
two jet rate compared to light quarks. For three jets, the situation is slightly different; 
for small value of y it is the four-jet rate that starts dominating. In this region therefore 
the three-jet rate for massive quarks starts to exceed the massless one, which tends to 
have a larger fraction of four-jet events. 

In Fig. [7J the results of the parton shower as implemented in APACIC++ are confronted 
with these analytical jetrates, which basically represent a resummed tree-level calcula¬ 
tion. The agreement over a wide range of phase space, down to small values oiy < 0.004 
is very good. Only for very small values of y the three-jet rate as obtained from the 
parton shower starts to overshoot the analytical result, which is due to the increased 
importance of four-jet events in this region. Finally, Fig. |H1 exhibits the ratio of three- 
jet rates with 6-quarks and with light quarks (m, d, and s) only. There, data from the 
Delphi collaboration IS2] are contrasted with a full next-to leading order calculation ESI 
with varying 6 quark mass, with the pure parton shower result as obtained by APACIC++ 
and with the result of the combined matrix element plus Sudakov weight. The data 
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Figure 4: Cross section for the production of up to four jets (left) and the correspond¬ 
ing Sudakov weight that will be attached in the merging procedure (right) vs. the jet 
resolution y. Two/three/four jet topologies are depicted in red/green/blue, massless 
jet-conhgurations are shown with solid lines, jet conhgurations involving two massive b 
quarks with rrib = 4.5 GeV are shown with dashed lines. 



iog,o(y) iog,o(y) 


Figure 5: Cross section times Sudakov weight vs. the jet resolution y for each jet 
conhguration and the sum of them. The left plot shows massless jets, the right one 
displays jet conhgurations with two massive b quarks. 
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Figure 6: Jetrates for two- and three jets from the combined matrix element plus 
Sudakov weight. 


clearly prefer a lighter b quark with 3 GeV in the NLO calculation ®. In addi¬ 

tion, the parton shower obviously describes the data surprisingly well; for the merged 
matrix element plus Sudakov weight result, a running b quark mass seems to be also 
appropriate. 


®This result is consistent with the value obtained by evolving the measurement at the T resonance 
up to the Z boson mass, cf. m 
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Figure 7: Analytical results for the jetrates from Fig. El vs. the results obtained from 
the parton shower. 



Figure 8: Heavy-to-light three-jet rate ratio i ?3 = in dependence on the 

Durham jet resolution y. Data (black points) from a Delphi measurement j221 are 
shown together with the SHERPA prediction (the dark green histogram corresponds to 
the shower result and the light green curve is obtained by combining matrix elements 
with Sudakov weights), and with an analytic calculation IHHI (blue lines, dashed=LO, 
solid=NLO) for two different values of mi, =3, 5 GeV. 
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Comparison of shower and hadron level: Hadronisation correc¬ 
tions 

Before further comparing results of the parton shower with experimental data, it must be 
stressed that many of such comparisons are feasible and meaningful on the hadron level 
only. This implies that the partons emerging from the shower must be hadronised with 
some phenomenological model. In the case of SHERPA, the choice at the moment is to 
employ the Lund string model for it. This model has a number of parameters to be tuned 
to data. The tuning of Monte Carlo event generators to data is an intricate procedure, 
involving an optimisation in a multi-dimensional phase space of parameters. These 
parameters may be perturbative (like, for instance, a^), or non-perturbative (like, for 
example, the string tension in the Lund model), or they could characterise the transition 
between the perturbative and the non-perturbative regime (e.g. the parton shower cut¬ 
off). For further details on such a tuning procedure, the reader is referred to j541 
The parton shower of APACIC++ in its versions 1.0 and 2.0 together with the Lund string 
model implemented in Pythia has been tuned in and E7|, respectively. 
Hadronisation and the tuning of phenomenological models, however, induces some source 
of systematic error in any simulation that needs to be investigated. In Fig. Elthe parton 
shower results for total jet rates in the Durham scheme for 2-5 jets are confronted 
with results after hadronisation. Evidently, both results coincide on a level of 10% or 
better down to jet resolutions of y ~ 0.001. This corresponds to relative transverse 
momenta of the order of 3 GeV, a kinematical regime, where decays of 6-hadrons start 
to matter. This also implies that down to such low values parton shower results can be 
compared directly to data. In Fig. HUl two event shape variables, thrust and oblateness, 
are studied. Again, results before and after hadronisation are compared. Pictorially 
speaking, thrust is a measure for how “jetty” an event is. A thrust of 1, for example, 
describes a perfect, “pencil-like” two-jet event, where both jets are oriented back-to-back 
and have no transverse spread. Clearly, such a conhguration is severely suppressed at 
the hadron level, since it corresponds to events with few, highly collinear particles only. 
Thus, in the T = 1 bin, hadronisation corrections are large. Over a wide range of thrust 
values away from such extreme conhguration, however, the shapes of the results before 
and after hadronisation are nearly identical; hadronisation corrections in this region 
therefore a just a constant that can be taken directly from data. In this region of, say, 
0.05 <T< 0.6, the parton shower results can be confronted directly with data. Similar 
reasoning holds true not only for the other exemplary observable, oblateness, but for a 
plethora of observables. 
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Figure 9: Durham jet rates at LEP 1. The shower level result (solid lines) of SHERPA is 
contrasted with its result after hadronisation (dashed lines). 



Figure 10: Thrust and Oblateness at LEP 1. The shower level result (solid lines) of 
SHERPA is contrasted with its result after hadronisation (dashed lines). 
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Comparison with experimental data from LEP 

Having validated the numerical correctness of the parton shower implementation, the 
focus shifts from the investigation of the behaviour of the Sudakov form factors to the 
study of experimental observables. There, experimental data can be confronted with 
the results obtained with the parton shower, eventually after a merging with multijet 
matrix elements. As mention above, in the previous section, the considered data are on 
the hadron level, consequently a hadronisation model has to be applied. The parameters 
of the model were tuned, the perturbative input consisted of a sample of multijet matrix 
elements for up to hve jets, merged with the parton shower of APACIC++. The quality of 
the tuning can be judged by considering the multiplicity distribution of charged hadrons 
and of their scaled momentum, cf. Fig. In this plot as well as in each of the following 
ones, data are confronted with the parton shower implemented in APACIC++(solid lines) 
and with the merged multijet matrix element plus parton shower (dashed lines), both 
after hadronisation. For the merged results, contributions from different jet multiplic¬ 
ities are indicated in different colours. Clearly, the merged sample tends to produce a 
slightly larger fraction of events with harder jets, leading to higher parton and hadron 
multiplicities; this is visible from the fact that the “parton shower only” sample slightly 
undershoots the bins with comparably high hadron multiplicity. For the momentum 
distribution, however, this minor trend washes out. In both cases, the results obtained 
from the simulation show excellent agreement with data, this is even more so when tak¬ 
ing into account that in all plots shown the statistical errors on the Monte Carlo results 
are of the same order as the experimental errors, indicated by the yellow error bands. 
The hrst real test of the parton shower performance is to check whether it is able to 
reproduce event shape observables, such as thrust, major, the C-parameter, or oblateness 
cf. Figs. 1121 and M In all cases, the agreement of the data with the generated events 
is excellent. In the “parton shower only” sample the trend mentioned above, namely of 
being a little bit softer than the merged sample, is continued. 

Turning from the event shape observables to jet observables, the small differences be¬ 
tween the two samples vanish nearly completely. In Fig. the relative fractions of 
events with different numbers of jets are exhibited in dependence on the jet resolution 
parameter y. The mutual agreement of the two samples and their agreement with data is 
excellent: The parton shower is perfectly capable of describing jet multiplicities at LEP 
1. This Ending is repeated for the diEerential jet rates, viz Fig. In Fig.^J the topo¬ 
logical structure of four jet events is investigated. To this end, the Bengtsson-Zerwas 
angle inn and the modified Nachtmann-Reiter angle |n21 are studied. Apparently, the 
merged sample is in perfect agreement with the data, whereas the “parton shower only” 
sample exhibits a slight shift away from them. This, however, is not a big surprise; after 
all, these observables do depend on interferences between different diagrams. To take 
this into account clearly is well beyond the abilities of the parton shower. 

^°For definitions, the reader is referred to appendix IQ 
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Figure 11: Charged multiplicity and scaled momentum at LEP 1. The left plot shows 
the distribution of the number of charged particles together with a measurement by 
Opal jSH]- The mean charge multiplicity is also stated together with its PDG value j59j . 
On the right hand side, a scaled momentum distribution is plotted against Delphi data 
M- The data are contrasted with results obtained through the parton shower alone 
(solid lines) with those obtained when the merging of matrix elements for up to hve jets 
and the parton shower has been employed (dashed lines). In the latter case coloured 
lines indicate the contributions from individual matrix elements. 
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Figure 12: Thrust and Major at LEP 1. The hadron level result of SHERPA is contrasted 
with measurement from the Delphi collaboration |M|. Line styles and colours are the 
same as in Fig. 
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Figure 13: C-Parameter and Oblateness at LEP 1. This plot shows the event shape 
variables C-parameter and oblateness, together with Delphi datajS^. Line styles and 
colours are the same as in Fig. El 
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Figure 14: Durham jet rates at LEP 1, taken from [Hn]- Line styles and colours are the 
same as in Fig. 
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Figure 15: Differential jet rates in the Durham scheme at LEP 1, taken from EDI. Line 
styles and colours are the same as in Fig. ITTl 



Figure 16: Four jet angle distributions. Shown are the Bengtsson-Zerwas angle (left) 
and the modified Nachtmann-Reiter angle (right). The data points are from a DELPHI 
measurement EH- 
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Tevatron 


Having investigated the final state parton shower, the focns shifts now on a stndy of the 
parton shower in the initial state. There, a good observable to judge the performance 
of the parton shower is the transverse momentum distribution of lepton pairs in Drell— 
Yan scattering. In Fig. El the p± distribution of such pairs with masses at the Z- 
pole (91 ± 15 GeV) is displayed. In order to describe the left side of the distribution, 
the initial partons have to be supplemented with intrinsic transverse momentum w.r.t. 
the hadron they stem from. In SHERPA, this intrinsic k± is distributed according to a 
Gaussian, with expectation value k± = 0.8 GeV. In this plot, two different simulations 
runs are confronted with data, both only use the parton shower for parton radiation. 
They differ in the choice of the starting scale of the shower evolution; one is starting 
at f = Mfi ^ M|, the other one is starting at the c.m. energy squared of the incoming 
hadrons. Glearly, this leads to differences in the treatment of hard radiation, and as 
expected in the former case, parton radiation ceases to exist at scales of the order of 
the hard scale. The agreement of both simulation runs with data at scales up to 40-50 
GeV, however, is excellent. In Fig. El the same observable is depicted once more, this 
time, however, the data are confronted with the pure parton shower, starting at the high 
scale, and with the merged result, including up matrix elements for the production of 
the Z accompanied with up to three hard jets. The results again are in great agreement 
with data. This shows that also the implementation of the initial state parton shower 
in APACIC++ has succeeded. 
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Figure 17: The p± distribution of the Z-boson in comparison with data from CDF at 
the Tevatron, Run I [^2] • The bottom plot shows the same distribution as the top one, 
but with focus on the low momentum region. The solid line indicates the shower result 
when using a fixed start scale of 1800 GeV, while the the dashed line is obtained when 
using s' ~ M| as start scale. The SHERPA results have been multiplied by a constant 
TC-factor of 1.45 to match the data. 
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Figure 18: The p± distribution of the Z-boson in comparison with data from CDF at 
the Tevatron, Run I |n21- This plot shows the same distributions as Fig. El but this 
time the pure shower performance is compared with the result obtained when merging 
matrix elements with up to 3 extra jets. The jet scale Qcut was hxed to 20 GeV. The 
coloured lines give the contributions from individual matrix elements. 
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5 Conclusions 


In this paper, the new version of the program APACIC++ has been presented. In contrast 
to its first version, APACIC++ ceased to be a stand-alone program; now, it is just the 
parton shower module of a larger framework. This transformation resulted in a reloca¬ 
tion of many steering and service classes. On the other hand, the scope of APACIC++ 
widened in such a way that in its present state it is also capable of performing the par- 
ton shower in the initial state of hadronic collisions. In addition, the merging procedure 
with the matrix elements has been refined. On the technical side, some algorithms have 
been changed, in particular, the look-up tables for the Sudakov form factors have been 
replaced by an algorithm based on the hit-or-miss method. 

Some exemplary results obtained with the new version of APACIC++ have been pre¬ 
sented as well. The comparison with known analytical results, that properly resum 
large logarithms, validates the physical and numerical correctness of the implementa¬ 
tion of the Sudakov form factors. This is not too obvious, since the analytical results 
contains logarithms in terms of transverse momentum scales, whereas the parton shower 
in APACIC++ has been formulated in terms of virtual mass, supplemented with appropri¬ 
ate scale choices and an explicit angular ordering. Both the parton shower alone and it 
being merged with matrix elements results in an excellent agreement with a collection 
of precise data. These results span up to five orders of magnitude and the agreement 
is, especially for the LEP 1 data, on the level of a few percent. This proves the validity 
and underlines the abilities of the parton shower implemented in APACIC++. 
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A Brief program documentation 

APACIC++ is the module responsible for parton showers inside the SHEEP A framework 
Therefore, basic physics tools, like e.g. four-momentum, PDFs, particle dehnitions, and 
the jet algorithm, are provided through the overall framework. Furthermore, a num¬ 
ber of features related to the jet-veto are closely connected to the implementation of 
the merging procedure in SHERPA. Nevertheless, APACIC++ could be employed by the 
framework of a different event generator, provided the basic physics tools are made 
available. 

This section gives a brief summary of the tasks each class in APACIC++ is responsible 
for. Where needed, some details on specihc implementation issues are presented that 
should, in principle, enable the interested user to implement and test some of his or 
her own ideas. As stated above, in Sec. 121 APACIC++ represents the evolving parton 
shower in terms of binary Trees, consisting of doubly linked Knots. Each of them 
represents one individual parton splitting. APACIC++ steers the shower evolution in 
the initial and hnal through two different classes, Initial_State_Shower and Final_- 
State_Shower, respectively. These shower classes hll the Knots through corresponding 
Sudakov form factors, encoded in Timelike_Sudakov and Spacelike_Sudakov, where 
the latter carries a link to the appropriate PDF. Therefore, there are two instances of 
Spacelike_Sudakov with potentially different PDF in the Initial_State_Shower. For 
convenience, both Sudakov classes are derived from the class Splitting_Group; as such, 
they contain all relevant Splitting_Functions. However, both showers reconstruct the 
branching kinematics of each splitting from quantities like the scale t and the energy 
splitting parameter This is achieved in the two classes Timelike_Kinematics and 
Spacelike_Kinematics, respectively. 

The interface with SHERPA: the class Apacic 

This class dehnes the general interface to the shower package APACIC++. The parton 
shower evolution for a set of given partons is performed along the following steps: 

• The parton shower of each event is started after setting appropriate initial condi¬ 
tions, taken from the merging with the matrix element. These initial conditions 

- basically starting scales t and maximal angles for the coherent shower evolution 

- are brought directly into the Trees by hlling the pseudo parton shower history 
of the matrix elements into corresponding Knots. Pointers to the Trees needed 
for this operation are extracted from the showers and handed to the outside world 
through the methods FinTreeO and IniTreesO. 

• By calling PerformShowersO, the corresponding methods in the class Final_- 
State_Shower and in the class Initial_State_Shower are triggered to perform 
the shower evolution. All necessary boosts at the beginning and at the end of the 
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shower^^ are carried out in here. The methods SetJetvetoPt2() and SetFacto- 
risationScaleO are used to pass settings relevant for the parton shower piece of 
the merging into the showers. The check of whether a jet was lost during shower 
evolution is performed through Final_State_Shower: :ExtraJetCheck(). 

• The hnal result of the shower evolution can be transfered to the SHERPA frame 
work with the help of ExtractPartonsO. 

Running the showers 

The class Final_State_Shower 

This is the central class of the hnal state shower implementation. It controls the se¬ 
quence of evolution steps^^. For the individual shower evolution of a single parton the 
class Timelike_Sudakov is employed. Branching kinematics are constructed after the 
evolution has hnished, using the class Timelike_Kinematics. It is also responsible for 
some kinematics checks during the shower evolution and for the jet veto. 

The hnal state shower evolution of a given jet ensemble is performed by Perf ormShower (), 
while the method FirstTimelikeFromSpacelikeO is called for the hnal state (time- 
like) shower of a parton emitted during the initial state (space-like) shower evolution. 
The diherent methods of the class Final_State_Shower are responsible for the following 
tasks: 


• The method Perf ormShower () initiates the hnal state shower on a given tree, i.e. 
starting from its root Knot. This is done in InitializeJetsO, the kinematics 
are constructed afterwards utilising TimelikeKinematics: :DoKinematics(). 

• InitializeJetsO initialises the jet system emerging from a given (dummy) 
mother knot. The algorithm is performed recursively along the following steps: 

1. If one or both of the two daughters are allowed to decay, the mother Knot 
is hlled accordingly and the parton system produced by the decaying daugh- 
ter(s) must be further evolved. These tasks are achieved by FillBranchO 
and EvolveJetO, respectively. 

2. If any the two daughters can not decay (internal lines of a ME), Initialize¬ 
JetsO is called again, with the daughter in question taking the role of the 
mother knot. 

• To £11 a branch, its splitting scale t and energy splitting variable have to be 
determined. In APACIC++, this is realised through FillBranchO. There, a mother 
with parameters {fa, Za} for a decay a ^ be into two already specified massless 
daughters serves as input for the determination of their {U, Zi} {i = b, c). 

^^For instance, the final state shower is always performed in the rest frame of the particles starting 
it. 

^^Note that the class Final_State_Shower is also utilised for the (time-like) evolution of any parton 
emitted off the initial state shower. 
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Step by step one of the daughters is chosen to obtain a new trial virtuality ti 
according to Eq. realised by Timelike_Sudakov: :Dice(). In each step, the 
last ti serves as starting condition for the determination of the new one. Only if 
the system of both daughters passes the kinematic constraints, this sequence of 
alternating reduction of the ti is interrupted and they are accepted. Then, the 
mothers energy splitting variable Za must be modihed according to Eqs. (naEi 
implemented in Timelike_Kinematics: : Shuffle() to compensate for the gain in 
virtuality of the daughters. 

• Evolve Jet 0 evolves a given parton system, consisting of a mother knot with 
two daughters. The algorithm works recursively: Both daughters may branch 
further, therefore FillBranchO is called first to determine both virtualities within 
the given kinematical constraints. In case neither of the daughters branches, the 
algorithm stops at once. Otherwise, EvolveJetO is called for the corresponding 
daughter(s). 

• FirstTimelikeFromSpacelikeO in contrast tries to initialise a new jet system 
emerging from a time-like particle emitted by a space-like shower. If, starting from 
the virtuality ta of the space-like branch b ac, the Timelike_Sudakov: :Dice() 
yields a suitable virtuality tc > to respecting the constraint Eq. (OH|l . daughters 
are initialised and the jet is evolved by EvolveJetO. 

• After the shower is completed, colours are set through the method SetAllColours () 
according to the algorithm described above. Then, the method ExtractPartons () 
extracts the partons (instances of the class Particle of the full framework) from 
the outgoing Knots of the Tree and hlls them into the event record. 

The class Initial_State_Shower 

The class Initial_State_Shower governs the space-like shower evolution. Starting with 
the particles entering the hard 2^2 piece of the process, a backward shower is per¬ 
formed with the help of the classes Spacelike_Sudakov (two objects, one for each beam) 
and Spacelike_Kinematics. For the treatment of time-like parton emissions the con¬ 
trol is transfered to the class Final_State_Shower. The lnitial_State_Shower object 
contains the following methods: 

• Perf ormShower 0 is called in order to start the space-like shower evolution. Infor¬ 
mation one the matrix element kinematics is given in form of two partially hlled 
trees, corresponding to the evolution of the left and the right incoming parton, 
respectively. Similar to the final state case, the shower is initiated by employing 
InitializeSystemO. After the evolution has finished some consistency checks 
are done, and the shower history is transfered back into the laboratory frame. 

• InitializeSystemO determines the initial system for the shower evolution. The 
algorithm starts with calling FillBranchO for each of the given matrix element 
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partons in order to obtain the virtuality t* and energy fraction Zi of the first split¬ 
tings. Then, the off-shell momenta are constructed with the help of Spacelike_- 
Kinematics: : InitKinematicsO. 

• Additional branchings are appended by EvolveSystemO. This methods performs 
the evolution of a given system of two space-like partons, by recursively appending 
additional branchings. In each step a system consisting two partons of different 
incoming beams. The shower evolution is performed backward, i.e. from the hard 
interaction towards the beam particle. The parton with larger virtuality is sup¬ 
posed to be closer to the hard interaction and is selected for the next evolution 
step. There, the virtuality t and energy fraction is hxed by consecutive calls 
of FillBranchO, CalculateMaxTO, and FirstTimelikeFromSpacelikeO. Af¬ 
ter each step, the four-momenta of the participating partons are evaluated with 
DoKinematicsO, leading to a new system of two space-like partons. 

• After the shower evolution is hnished, the created parton set can be accessed by 
calling ExtractPartonsO. It hlls all hnal state particles connected to the two 
incoming Trees into the event record. 

• Two conhguration methods need to be mentioned, which are of special importance 

in the context of merging ME and PS, namely Set JetvetoPt2() and SetFacto- 
risationScaleO. In order to guarantee a clean separation between, the jet-veto 
is imposed on any trial emission inside the shower. Usually, the p± for any trial 
radiation is restricted by the jet resolution scale For the matrix element 

with the highest multiplicity of jets, the separation cut is dynamically hxed by the 
smallest transverse momentum present in the hard interaction, owing to the fact 
that shower is not supposed to produce radiation harder than any QCD radiation 
present in the matrix element. In full analogy, merging ME and PS involves 
evaluating the PDFs at a specihc factorisation scale pp- The dependence on the 
Pf is to be cancelled by a corresponding PDF evaluation during the hrst parton 
emission in the initial state shower. The factorisation scale pp does not necessary 
coincide with the jet-veto scale mentioned above p™*. For instance for (leptonically 
decaying) W bosons produced in hadronic collisions, there is no jet-veto applied, 
but of course pp is not vanishing. In fact, in this example case, the factorisation 
scale is identihed with the hard c.m energy . 

Splitting functions &; Sudakov form factors 

The class Splitting_Function is purely virtual. It dehnes a common interface to the 
splitting functions Pa^bc{.z). Among others, this class and its specihc instantiations 
include methods to 

• access its havours (GetFlAO, GetFlBO, and GetFlCO) 

• to determine according to an approximative distribution (GezTO), and its inte¬ 
gral in a given region (CrudeIntO), 
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Figure 19: class hierarchy 


• a correction weight (GetWeightO), and the exact functional form of the splitting 
function through its operator()) 

• and some methods to manage Splitting_Groups. 

All implemented splitting functions are derived from this class, cf. Tab. 1201 This 
abstraction renders the inclusion of further splitting functions an easy task. 

The class Splitting_Group is a container derived from the abstract class Splitting_- 
Function. It contains all splitting functions of a given flavour, either in the forward 
splitting or for backward splitting as a Backward_Splitting_Group object. It is thus 
responsible for the determination of a branching for one specific flavour. The integration 
routine therefore returns the sum of the integrals of all single Splitting_Functions. 
After the integration the routine SelectOneO can be called to choose one splitting 
mode out of the available options according to the integrals. All subsequent calls to 
flavour access methods, dice routines or the weight calculation then correspond to the 
selected branching. Note that the Splitting_Group is also the base of any Sudakov 
form factor determination, as depicted in Fig. 

Both Sudakov classes are derived from the corresponding Splitting_Group. Of course, 
in the space-like case, the mother class is the Backward_Splitting_Group. Following 
their names, the classes Timelike_Sudakov and Spacelike_Sudakov govern the time- 
like and space-like shower evolution, respectively. The latter one must be supplemented 
with a link to a PDF implementation, realised in SHERPA through an object derived from 
an abstract PDF_Base. However, both Sudakov classes store appropriate Splitting_- 
Groups for every flavour taking part in the shower. In both classes, the method DiceO 
is responsible for the determination of virtualities (Produce! ()), daughter flavours, and 
energy fractions. The algorithm chosen is the hit-or-miss method, and both classes 
thus contain various veto methods implementing cuts and correction weights. To exem¬ 
plify this, both classes have a method CplVetoO incorporating the correction weight 
cXsiQ)/(^s{.Qmin)■ Especially for the incorporation of issues related with coupling con¬ 
stants, both have a link to the class Sudakov_Tools. 
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QED splitting functions 



Figure 20: A summary of implemented splitting functions and their class names in 
SHERPA. 

Kinematics 

Kinematics are implemented in two classes, Timelike_Kinematics and Spacelike_- 
Kinematics for the determination of the kinematics in the hnal and initial state showers, 
respectively. They include 

• checks whether a branching is kinematically possible (KinCheckO), 

• a check whether an emission yields a jet that must be vetoed in the merging 
procedure (JetVetoO), and 

• a method to construct kinematics (DoKinematicsO). 

In addition, Timelike_Kinematics provides a method to redefine the energy splitting 
according to Eqs. (HUES) (ShuffleO), and Spacelike_Kinematics has a method 
to determine the maximal kinematically allowed t for the initialisation of a final state 
shower off an initial state splitting (CalculateMaxTO), cf. Eq. (0!^ 

Basic structures 

Trees are used as representations of parton shower histories in terms of interconnected 
binary splittings. It provides the basic structure and all necessary routines to handle 
operations, like boosts, on the whole or parts of a tree. Three trees form a complete 
parton shower history: Two trees correspond to the initial state shower evolution of the 
right and the left incoming particle, respectively, and the other tree holds information 
of the final state shower of the outgoing particles, cf. Sec. El 

The Knots are the basic elements forming a binary Tree. They store all properties of a 
single branching, like flavour, momentum, virtuality, energy component, etc.. 
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B Altarelli—Parisi splitting functions 

In this appendix all splitting fnnctions, relevant for the shower evolntion of APACIC++ 
are listed, cf. jnn insi- The splitting fnnctions are obtained after averaging over the 
azimuthal angle. The mass terms in each case are parametrised by the variable = 
{nil + rn^)/[{pi + PjY — ''^(p)]- The expressions for splittings with the decay products 
exchanged fulhl the obvious symmetry relation Pa^bc{,z) = Pa^chi)- — z). 

Fermions and photons: 



(33) 


(34) 


Quarks and gluons: 



(35) 

(36) 

(37) 


7)7^33(2) — Za ^-h 2(1 - 2) 

1 . z z 


In the massless limit (/i —0), these splitting functions reduce to the well-known form 


1 + 

P„(^) = Ca , 

P„( 2 ) = T„ [2" + (1 - 2)"] , 


(38) 

(39) 

(40) 



Gluinos and gluons: 



(41) 

(42) 


Pg^gg(z; /i^) = Ca [z^ + (1 " zf + . 


Squarks and gluons: 



(43) 


(44) 
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C Definitions of event shapes 

The global properties of hadronic events may be characterised by set of observables, 
usually called event shapes. In section HJ the following shape observables have been 
considered. 


• Thrust T : 

The thrust axis Ht maximises the following quantity 


T 


max 

nx 


f Ei \Pi-nT\ \ 

V EM\ 


(45) 


where the sum extends over all particles in the event. The thrust T tends to 1 
for events that has two thin back-to-back jets (“pencil-like” event), and it tends 
towards 1/2 for perfectly isotropic events. 


• Thrust Major TMajor : 

The thrust major vector nMajor is dehned in the same way as the thrust vector, 
but with the additional condition that UMajor must he in the plane perpendicular 
to Ht- 


TMajor = _ max__ 

‘^M.ajor EtIx 


Eli \P'i ' ^Major | 

Ei \Pi\ 


(46) 


• Thrust Minor Tninor : 

The minor axis is perpendicular to both the thrust axis and the major axis, 
^Minor = ut X UMajor- The value of thrust minor is then given by 


TMinor 


Eli \Pi ' ^Minor 

Ei IS I 


(47) 


• Oblateness O : 

The oblateness is defined as the difference between thrust major TMajor and thrust 
minor TMinor : 

O = TMajor — TMinor (48) 


• C-parameter C : 

The C-parameter is derived from the eigenvalues of the linearised momentum 
tensor 0 “^ , defined by 


0«/3 


1 PfPi 

Ei \Pi\ i \Pi\ 


a, 13= {x,y,z} . 


(49) 


The three eigenvalues A* of this tensor dehne C with 


C — 3 (A 1 A 2 + A 2 A 3 -|- A 1 A 3 ). 


(50) 
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